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ABSTRACT 

Narrow depolarized canals are common in maps of the polarized synchrotron emission of 
the Milky Way. Two physical effects that can produce these canals have been identified: the 
presence of Faraday rotation measure (RM) gradients in a foreground screen and the cumu- 
lative cancellation of polarization known as differential Faraday rotation. We show that the 
behaviour of the Stokes parameters Q and U in the vicinity of a canal can be used to identify 
its origin. In the case of canals produced by a Faraday screen we demonstrate that, if the polar- 
ization angle changes by 90° across the canal, as is observed in all fields to-date, the gradients 
in RM must be discontinuous. Shocks are an obvious source of such discontinuities and we 
derive a relation of the expected mean separation of canals to the abundance and Mach num- 
ber of supernova driven shocks, and compare this with recent observations bv lHaverkorn et alJ 
(2003). We also predict the existence of less common canals with polarization angle changes 
other than 90°. Differential Faraday rotation can produce canals in a uniform magneto-ionic 
medium, but as the emitting layer becomes less uniform the canals will disappear. We show 
that for moderate differences in emissivity in a two-layer medium, of up to 1/2, and for Fara- 
day depth fluctuations of standard deviation < 1 rad, canals produced by differential rotation 
will still be visible. 
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1 INTRODUCTION 

Polarized (synchrotron) radio emission is a rich source of in- 
formation about the relativistic and thermal plasmas and mag- 
netic fields in the interstellar medium (ISM). Recent observa- 
tions have revealed an abundance of unexpected features that 
arise fro m the propagation of the emission through the turbu- 
lent I SM I Wier inga et alll993l:IUvaniker et alll998l:lDuncan et alJ 
| !999l:lGrav et alll999l:lHaverkorn et all200d:lGaensler et alJ20oTF 
Woll eben et alJ l200rj) . Arguably the most eye-catching is the 
pattern of depolarized canals: a random network of dark nar- 
row regions, clearly visible against a bright polarized back- 
ground. These canals evidently carry information about the ISM, 
but it is still not quite clear how this information can be ex- 
tracted. Two theories for the origin of the canals have been pro- 
posed; both attribute the canals to the effects of Faraday ro- 
tation, but one invokes steep gradients of the Faraday rotation 
measure (RM) across the telescope beam, in a Faraday screen 
jHaverkorn. Katgert & de BruvnlEoOcL Eo04). whereas the other 
relies on the line-of-sight effects pro ducing differential Faraday ro- 
tation (Beck 1999; Shukurov & Berkhuiisen 2003). In order to use 
the canal properties to derive parameters of the ISM, one must cor- 
rectly identify their origin. 

In this letter we briefly discuss a few important aspects of 
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the two theories describing the origin of canals. Detailed dis- 
cussion of relevant depolar ization mechanisms is presented in 
iFletcher & ShukuroU l2006aL and in preparation). Canals produced 
by Faraday rotation measure gradients are discussed in Section |2| 
where we show that these canals require a discontinuous distribu- 
tion of free electrons and/or magnetic field; we further suggest an 
interpretation of the discontinuities in terms of interstellar shocks. 
The case of differential Faraday rotation is discussed in Section^ 
in Section |4| we suggest an observational test that can be used to 
identify the specific mechanism that produces a given canal. 
The defining features of the canals are as follows: 

(i) the observed polarized emission P approaches the polariza- 
tion noise level ap, P < ap; 

(ii) the canal is about one beam wide; 

(iii) the canal passes through a region of significant polarized 
intensity, say P > 3a p; 

(iv) the canal is not related to any structure in total intensity, 
and so cannot be readily explained by, e.g., an intervening gas or 
magnetic filament. 

Polarized intensity vanishes when emission within the telescope 
beam consists of two equal parts with mutually orthogonal polar- 
ization planes. Thus, feature [(I)| most often arises because the po- 
larization angle \£ changes by 90° across the canal. However, in 
Section|2|we predict a new type of canal across which the observed 
polarization angle does not change. 
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We only consider cana ls occurring in properly calibrated 
maps: Haverk orn et alJ <Eo04) argue that the observations we dis- 
cuss in Section lTTl do not suffer from missing large-scale structure; 
Reich ( 2006) discusses the calibration of radio polarization maps in 
depth. 

Polarized radiation is commonly described in terms of the 
complex polarization, 



F i 



V = pexp(2i^) 



(1) 



where p, the degree of polarization, is the fraction of the radia- 
tion flux that is polarized. When polarized emission passes through 
magnetized and ionised regions, the local polarization angle ip (at 
position r) changes by an amount depending on the wavelength A 
due to the Faraday effect 



ip(r) = ipo(r) + < 
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where K = 0.81rad m -2 cm 3 /iG -1 pc _1 is a constant, n e is the 
number density of free thermal electrons, B z is the component of 
the magnetic field along the line of sight (here aligned with the z- 
axis), and the observer is located at z — > oo. <j>{z) is known as 
the Faraday depth to a position z and gives the change in polariza- 
tion angle of a photon of wavelength A as it propagates from z to 
the observer. The maximum amount of Faraday rotation in a given 
direction is called the Faraday depth 1 
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The observed amount of Faraday rotation, determined by the ro- 
tation measure RM = d^/dA 2 , cannot exceed F, i.e., |RM| ^ 
|F|A" 2 . 

The value of RM is related to F, but often i n a complicated 
manner (see, e.g., Burn 1966; Sokoloff et al. 1998). Simplest is the 
case of a Faraday screen, where the source of synchrotron emission 
is located behind a magneto-ionic region (e.g., because relativistic 
and thermal electrons occupy disjoint regions): then RM = FA -2 . 
In a homogeneous region, where relativistic and thermal electrons 
are uniformly mixed, RM = 0.5FA -2 . 

Observations of linearly polarized emission provide the Stokes 
parameters I, Q, U which are related to p and \£ via V — (Q + 
iU)/I: 



{Q 2 + U 2 )^ 2 



U 



arctan — — 7^71" (sign Q — 1) sign U 



(2) 
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and the polarized intensity is P — (Q 2 + U 2 ) 1 ^ 2 = pL The com- 
plex polarization can be written in terms of 0, 



(4) 



V=^ [ W(r±)e(r)exp{2i[Mr) +t(r)]}dV , 
1 Jv 

where integration extends over the volume of the telescope beam 
V, W(r±) defines the shape of the beam, a function of position in 
the sky plane r± — (x,y), and e(r) is the synchrotron emissivity. 
The total intensity is similarly given by / = J y W(r±)e(r) dV\ 
the Faraday depth is a function of r±, F = F(r±). 



1 Thi s ter minology may cause confusion: many authors, including iBurnl 
i 19661) and lSokoloff et aljJl998L define the Faraday depth as F/X 2 , in our 
notation. However, it is more convenient, and physically better motivated, to 
define the Faraday depth, sim ilarly to th e optical de pth, as a dimensionless 
quantity, as used bv lSpangleJ Il982h and lEilekl j 1989 ) 




Figure 1. Sketch showing the variation in Faraday depth F with respect to 
transverse distance in the sky plane x. The Faraday depth changes by AF 
across the distance Ax and by AFd within the beam width D. 



2 CANALS PRODUCED BY A FARADAY SCREEN 

Consider polarized emission from a uniform background source 
passing through a Faraday screen - a layer where no further emis- 
sion occurs, but which rotates the polarized plane. If F varies with 
rj_, adjacent lines of sight within the beam are subject to different 
amounts of Faraday rotation and the observed degree of polariza- 
tion decreases. If the variation of F within the beam, AFd (Fig.[T), 
produces a 90° difference in ^, one might expect that the depolar- 
ization will be complete as illustrated in Fig.|2ja), and a canal will 
be observed along contours defined by AFd = (n + l/2)ir with 
n = 0,l,2,.... 

The situation is more subtle, though. Figure |2[b) shows ^ 
changing smoothly by 90° across a beam, as a result of a mono- 
tonic change in Faraday depth by the same amount, AFd = tt/2, 
as in panel (a). It is obvious that polarized emission does not can- 
cel within the beam; the polarized intensity will be lower than for 
a uniform arrangement of angles but there will still be a polarized 
signal detected with a polarization angle of about —45°. 

Thus, gradients in F across the beam can produce complete 
depolarization in a Faraday screen if F changes by AFd — 
(n + 1/2) 7r within a small fraction of the beam width. As shown 
in Fig.|3](solid line), obtained using Eqs (|2j and (0, the depth of a 
canal will only be less than 10% of the surrounding polarized emis- 
sion if the gradient in F occurs over one fifth or less of the beam: 
Ax/D < 0.2, where Ax is the extent of the region where the 
gradient occurs and D is the beam width. lHaverkorn et alJ (2004) 
reached a similar conclusion in an analysis of the depth and profile 
of canals observed at A84 cm (see their Sect. 4.4). 

As illustrated in Figs.|3c) and [3] (dotted line), a variation in 
F of a magnitude AFd = tt across the FWHM of the Gaussian 
beam (i.e. AF = tt and Ax/D ~ 1) can also produce strong 
depolarization but now the change in polarization angle across the 
beam is A^ c± 0°. Similarly other gradients in F can significantly 
reduce the degree of polarization; for example AF = 0.7tv (dashed 
line in Fig. [3j across a region that is about 0.7 x the beam width 
will produce complete depolarization and a change in angle of 54° 
across the beam. 

Larger gradients in F also result in strong depolarization and 
when AF > n the degree of depolarization is less sensitive to 
the resolution. This is illustrated in Fig. [3] for the case where 
AF = 1.57T (dash-dotted line); complete depolarization occurs 
when Ax/D = 0, 1.5 but in the range < Ax/D < 2 the degree 
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Figure 2. An illustration of depolarization in a Faraday screen: the shaded circle indicates the telescope beam, and dashes show the orientation of the 
polarization plane at various positions, (a) Abrupt change of the polarization angle \£ by 90° in the middle of the telescope beam leads to complete cancellation 
of polarization. This can be caused by a discontinuous change in F by 90° in a Faraday screen, (b) A similar but continuous change of ^ does not result in 
strong depolarization, (c) A continuous change in ^ by 180° can produce a canal. 
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Figure 3. The minimum degree of polarization in a canal produced by a 
gradient of Faraday depth in a foreground screen, p m i n , as a function of the 
relative linear extent of the gradient Ax/D, where D is the FWHM of a 
Gaussian beam and x is the position in the sky-plane (see Fig.[l). Different 
curves represent different increments in Faraday depth, AF (identified in 
the legend). For Ax/D < 1 (an unresolved gradient), F changes over a 
distance smaller than the beam width and AF is effectively discontinuous 
for Ax/D < 1. 



of polarization remains below the 10% level. Thus an increment of 
AF = 1.57T, and more generally AF = (n + l/2)7r, will produce 
one beam wide canals when the gradient occurs in a region roughly 
equal to or less than the beam width, i.e. where Ax/D < 1. Other 
large gradients in F, such as AF = nir, depolarize to around 10% 
of the background p and therefore can also produce canals; how- 
ever in this case one beam wide canals will only be formed when 
Ax/D ~ 1. 

So, if F is a continuous, random, function of position (with 
sufficiently strong fluctuations on the scale of the beam) we pre- 
dict the existence of canals with every change in polarization angle 
0° < A^ < 90° across the beam. On the other hand a discon- 
tinuous change in F by (n + l/2)ir, will produce canals across 
which the polarization angle changes by A\£ = 90°. This insight 
is important as to-date all observed canals have A^ c± 90°. The 
only other way in which Faraday screens can consistently produce 
canals with this change in angle is if all the gradients in Faraday 
depth on the approximate scale of the beamwidth are of the mag- 
nitude AF ~ (n + l/2)7r, n = 1, 2, . . . (n = in this case only 
produces weak depolarization, as noted above); this is extremely 
unlikely in a random medium. 

If a system of observed canals has A^ = 90° and is caused 
by a foreground Faraday screen — see Section |4| for a diagnostic 
test for the origin of canals — then the distribution of rotation mea- 



sure in the screen mu s t be di scontinuous on the scale of the beam. 
Haverkor n & Heitschl fe004l) used MHD simulations of Mach 10 
turbulence to show that RM gradients can be steep enough to pro- 
duce depolarized canals on smoothing with a sufficiently large sim- 
ulated beam. Note that a finite beam size is essential for the produc- 
tion of canals by this mechanism; with perfect resolution the width 
of the shock front will be resolved and no canal will appear. Fur- 
thermore, these canals will not fill up if the data are smoothed (solid 
and dash-dotted curves for Ax /D — > in Fig. [3} whereas canals 
produced by gradients other than AF = (n+ 1/2)tt will disappear 
under smoothing at a rate shown with dashed and dotted curves in 
Fig.0 

2.1 The mean separation of canals produced by shocks 

An obvious source of Faraday depth discontinuities are shocks in 
the ISM. Then the mean distance between the canals produced 
by a Faraday screen will be related to the distance between the 
shock fronts. The expected separation of suitable shock fronts can 
be derived using the model of interstellar shock- wave turbulence 
of iBvkov & Tontvginl fl987h and then compared with the separa- 
tion of canals observed by Haverkorn, Katger t"& de Bruvr] (12003). 
claimed to arise in a Faraday screen. 

If the frequency distribution of interstellar shocks of Mach 
number A4 is G(M), the mean number of shocks of a strength 
exceeding M. that cross a given position in the ISM per unit time 
can be obtained as 



N(M) 



r oo 
J M 



G(M)dM 



(5) 



and the mean separation between the shocks in three dimensions 
follows as 

C s 



L(M) 



(6) 



N(M) ' 

where c s is the sound speed. The separation in the plane of the sky 
is reduced by 7r _1 J* sin 2 OdO = 1/2 due to projection effects and 
multiplied by a further factor L/d if the shocks occur in a region of 
depth d. At an average distance to the shocks of d/2, the average 
angular separation in the sky plane is then 



d 2 ' 



(7) 



For supernova-driven shocks, Bykov&Toptysin (1987b de- 



rived 

Q(M) = G [M 



-(a+l) 



+ 3C(a)f cl (M 



(8) 



allowing for both primary and secondary shocks, where a is a 
numerical factor (a = 2 for a supernova remnant in the Sedov 
phase, and a = 4.5 for a three-phase ISM); C(a) w 2.3 x 10~ 2 , 
4.1 x 10 -3 for a = 2, 4.5 respectively; f c \ is the volume filling fac- 
tor of diffuse clouds that reflect the primary shocks to produce the 
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secondary ones; Go = aS^irrQ with S the supernova rate per unit 
volume and ro the maximum radius of a primary shock. We then 
obtain the following expression for the separation of shock fronts 
with M ^ M* in three dimensions: 



L(M) 



)pc( 
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15 kpc 
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h \ ( r 



50 pc/ V 10 °P C 



_J aC(a)f c \ 

M% (M* - l) 3 



(9) 



where z^sn is the supernovae rate and R and h are the radius and 
scale height of the star-forming disc. The term in square brackets is 
approximately 1 for a = 4.5 and M* — 1.2. 

If the magnetic field is frozen into the gas and the gas density 
increases by a factor e p at the shock then so will the field strength 
(we neglect a factor ~ 1 / y/2 due to the relative alignment of the 
shock and field). A canal forms where F increases discontinuously 
by AFd = (n + l/2)7r. So the required density compression ratio 



+ 1, 



(10) 



where (F) is the mean value of F. The gas co mpression ratio de- 
pends on the shock Mach number M as (e.g. JLandau & Lif shitz 
1960) 



M 2 



( 7 +l)M 2 r 
(7- l)M 2 + 2 ~' ± M 2 + 3 



(11) 



where we have used 7 c± 5/3 for the ratio of specific heats, which 
yields the value of the shock Mach number required to produce a 
canal: 



M* 



3e p 



1/2 



-1/2 



v 3(|AiV(F>| + l)V2 
In the field of canals observed by lHaverkorn et alJ j2l 



(12) 



at 



A84cm, (RM) c± — 3.4radm -2 and a by-eye estimate of their 
mean separation gives C c± 45'. The most abundant canals are 
produced by the shocks which can generate AFd = tv/2 and 
from Eq. IT2t this requires M* c± 1.2. Using / c i = 0.25, tak- 
ing a = 4.5 and the parameter values used to normalize Eq. J9}, 
we obtain L c± 10 pc for shocks with .M ^ .M*. Our estimate of 
L includes shocks with .M > .M* that will not produce canals, but 
the strong dependence of Q in Eq. J8} on Mach number means that 
L will underestimate the separation of canal generating shocks in- 
significantly. Using Eq. {71 , with L ~ 10 pc and C ~ 45', we find 
that the canals observed by Haverkorn et al. (2003) are compatible 
with a system of shocks occurring in a Faraday screen with a depth 
of d c± 100 pc. The maximum distance in this field, beyond which 
emission is comp letely depolarized, is estimated to be ~ 600 pc by 
Haverkor n et alJ fc003l) . Thus this model for the canals' origin re- 
quires that the nearest 100 pc in the direction (/, b) = (161°, 16°) 
is effectively devoid of cosmic ray electrons; if the cosmic rays nor- 
mally spread over a distance vxt ~ 1 kpc, where va — 20 km s _1 
is the Alfven speed and t c± 3 x 10 7 yr their lifetime, such a con- 
dition is difficult to explain. 



3 CANALS PRODUCED BY DIFFERENTIAL FARADAY 
ROTATION 

In Section[2]we discussed the depolarization effects of Faraday ro- 
tation measure gradients transverse to the line of sight acting on 
smooth polarized background emission. Now we will consider a 
uniform layer in which both emission and Faraday rotation occur. 
This g ives rise to t he well known effec t of differential Faraday ro- 
tation (Burn 1966; Sokolo ffetaill99Sh . where polarized emission 
from two positions along the line of sight separated by a rotation 
of A0 = 7r/2 exactly cancel, thus reducing the degree of polariza- 
tion. When a line of sight has a total rotation of F = mr there is 
total cancellation of all polarized emission in the layer and the de- 
gree of polarizat ion is p = 0. Canals produced b y this mechanism 
are discussed by Shukurov & Berkhu iiserl fc003l) . 



3.1 Differential Faraday rotation in a non-uniform medium 

We now investigate how canals produced by differential Faraday 
rotation are affected by deviations from uniformity in the magneto- 
ionic medium. First we will discuss the case of a two-layer medium 
in which the synchrotron emissivity is different in each layer. Then 
we allow for random fluctuations of F (the latter case produces 
what is known as Faraday dispersion). 

For a two-layer medium the fractional polarization can be 
written as 
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( 2 sin 2 Fi 2 sin 
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smFi smF 2 \ 
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(13) 



where we have assumed ipo = 0, h and 1 2 are the synchrotron 
fluxes from the furthest and nearest layers to the observer respec- 
tively, / = /1 + / 2 , and Fi and F 2 are the Fa raday depths throug h 
the two layers. [This is similar to Eq. (10) in lSokoloff et al J (11998). 
but with typos corrected.] We assume that F — F± + F2 = 7r, 
the condition for canals to form due to differential Faraday rota- 
tion, and investigate what will happen to the canals if h / I2. We 
parametrize the difference in the synchrotron emissivity in the two 
layers as e = (h — /i)//i, choose equal Faraday rotation through 
each layer for simplicity so F± — F2 — 7r/2, and substitute these 
into Eq. J13I . We then obtain the dependence of the minimum de- 
gree of polarization on e: 



1 



po 



7T (l + e/2) 2 



71 



(14) 



Thus, for a moderate difference in emissivity of e — 1/2 between 
two layers the canals will have a depth of p/po ~ 1/6 compared to 
p/po — for a strictly uniform medium; these canals will still be 
clearly visible and can be interp reted as e.g. contours of F or RM 
I Shu kurov & BerkhuiiserJl2003h . Where the difference in emissiv- 
ity is greater, say e ~ 2 or more as one might expect viewing distant 
bright spiral arm emission near the Galactic plane, the canals will 
more readily fill up, become much less distinct and more difficult 
to interpret confidently. 

Now let us consider the effect of random fluctuations in Fara- 
day rotation, so metimes called intern al Faraday dispersion. We start 
with Eq. (34) of lSokoloff etaHil998h : 



V = p 



l-exp(-S) 



(15) 



where S — 2a 2 ? — 2iF and gf is the standard deviation of the 
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Faraday depth F. In order to produce canals we need F — tt and 
then 

a F «l. (16) 

The term in brackets shows the powerful depolarizati on res ulting 
from internal Faraday dispersion fsee ISokoloff et alJ f° r 
large enough crp depolarization is complete. However, as long as 
cff < 1, canals produced by differential Faraday rotation will only 
fill up by 1/3 or less. For example, at A20 cm the canals will not 
be destroyed as long as the dispersion in Faraday rotation measure 
is less than 30radm~ 2 This is why canals are s till vi sible in the 
A20 cm map studied by Shukurov & Berkhuijsen ( 2003) which has 
a dispersion in rotation measure of about ctrm — 10 rad m -2 , i.e. 
(J f — crm A 2 ~ 0.4. 



(iii) If shocks produce the discontinuities in a foreground Fara- 
day screen that generate canals, the mean separation of the canals 
can provide information about the Mach number and separation of 
shocks in the screen fSection l2Jl . 

(iv) Canals produced by differential Faraday rotation are sensi- 
tive to non-uniformity in the medium along the line of sight, sys- 
tematic or random. However, they will remain recognisable if the 
synchrotron emissivity varies by less than a factor of about 2 or 
if the standard deviation of the Faraday depth is cff < 1 (Sec- 
tionlTTl. 
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4 OBSERVATIONAL DIAGNOSTICS FOR THE ORIGIN 

OF A CANAL 

The complex polarization emitted b y a layer p roducing differential 
Faraday rotation can be written as (Burn 1966) 



Po- 



sin F 



2i(^ + F/2) 



(17) 



We can see immediately that if ipo =0 we have Q/I = ReV oc 
sin F cos F and U/I = Im V oc sin 2 F. Since the orientation of 
the coordinate system in which Q and U are defined is arbitrary 
(i.e. the reference line from which we measure polarization angles 
can have any orientation), we can always choose a system in which 
ipo = near a canal. At the axis of a canal Q = U = 0, but in 
the case of a canal produced by differential Faraday rotation there 
exists a reference frame in the (Q, U) plane where one of the two 
Stokes parameters changes sign across a canal but the other does 
not. 

If, otherwise, a canal is produced by disc ontinuities in Faraday 
rotation across the beam (Sect. |3 we have (Fletcher & Shukurov 
2006a) 



'+ [U _ sinAi^ 2i(^ 0+ F) 
I _ AFd 



(18) 



in the vicinity of a canal and both Q and U will change sign across a 
canal (given ip = 0) except in the special case AFd = F + 2nix. 
The change in sign occurs since Q/I = pocos[2('0o + F)] on 
one side of the canal and Q/I = po cos [2(^0 + F) + tt] on the 
other: a similar variation occurs in U/I. Therefore the behaviour of 
the observed Stokes parameters Q and U across a canal provides a 
way to distinguish between canals produced by foreground Faraday 
screens and those resulting from differential Faraday rotation. 



5 SUMMARY 

The main points of this paper can be summarized as follows: 

(i) The behaviour of the Stokes parameters Q and U in the vicin- 
ity of a canal allows one to identify whether a foreground Faraday 
screen or differential Faraday rotation is the cause of the canal (Sec- 
tionE). 

(ii) A foreground Faraday screen can produce canals with any 
polarization angle change across the canal, < < 90°. How- 
ever, discontinuous jumps in the Faraday depth will only produce 
canals with A^ ~ 90° (Section|2j. 
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